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^ ! Abstract 

The evolution of electromagnetic and massless Dirac perturbations in the space-time geometry 



of Kehagias-Sfetsos(KS) black hole in the deformed Hofava-Lifshitz(HL) gravity is investigated 
and the associated quasinormal modes(QNMs) are evaluated from time domain integration data. 
We find a considerable deviation in the nature of field evolution in HL theory from that in the 
Schwarzschild space-time and QNM region extends over a longer time in HL theory before the 



5h . 

power-law tail decay begins. The dependence of the field evolution on the HL parameter a is also 
£Sj , studied. In the time domain picture we find that the length of QNM region increases with a. But 

0^ the late time decay of field follows the same power-law tail behavior as in the case of Schwarzschild 

■ black hole. 
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I. INTRODUCTION 



Recently a renormalizable theory of gravity in 3+1 dimensions was proposed by Hofava, 
inspired from the Lifshitz theory in solid state physics, now known as Hofava-Lifshitz(HL) 
theory[l-3]. The theory is a potential candidate of quantum field theory of gravity. It 
assumes a Lifshitz-like anisotropic scaling between space and time at short distances, char- 
acterized by a dynamical critical exponent z = 3 and thus breaking the Lorentz invariance. 
While in the IR limit it flows to z = 1, retrieving the Einstein's General Relativity(GR). The 
discussions on the consequences of theory are going on (4, 5]. As a new theory, it is interesting 
to investigate its various aspects in parallel. Since the theory has the same Newtonian and 
post Newtonian corrections as those of GR, systems of strong gravity, like black holes, are 
needed to get observable deviation from the standard GR. Various black hole solutions are 



found in HL theory 



25J of which one with asymptotically flat Minkowski spacetime is the 



KS black hole obtained by applying deformation in the original theory 
of KS black hole were explored in the past 



27 



261 ] . Various aspects 



39]. 



Perturbation of black hole is the best way to study the properties of these exotic objects. 
As it is well established, the evolution of perturbations in black hole space times involves 
three stages 4CH42I]. The first one is an initial response containing the information of the 
initial form of the original wave field followed by a region dominated by damped oscillation 
of the field called quasinormal modes, which depends entirely on the background black hole 
spacetimes. Finally the field decays at late time as a power-law fall off of the field. QNMs 
play an important role in astrophysical processes involving black holes. These interests are 



associated with their relevance in gravitational wave analysis [43|]. Apart from the obser- 
vational interest of detection of quasinormal ringing by gravitational wave detectors [44]. 
study of QNM of AdS black holes are of particular interests in the study of AdS/CFT 



correspondence 



45|, 



461 ] . It was conjectured that a large static black hole in AdS space corre- 



sponds to a thermal state in conformal field theory (CFT). The fundamental link between 
QNMs of a large AdS black hole and the dual field theory were shown in the work of Horowitz 
who computed the QNMs of AdS black holes 47j] . Motivated by this, a series of studies on the 



QNMs in asymptotically AdS backgrounds were done for various fields and dimensions 
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511 ] . See Ref. 52] for a recent review of QNMs of the black hole. The late-time relaxation 



of field perturbation is also an interesting topic of study since it reveals the actual physical 



2 



mechanism by which a perturbed black hole sheds its hairs. It was first demonstrated by 
Price that in the background of the Schwarzschild spacetime the massless neutral field at 
late time dies off as \l/ ~ t^( 2f + 2 ) or \I> ~ £-( 2 ^+ 3 ) depending on the initial conditions and the 



multipole order 



531 ] . The evolution of field and the spectra of QNMs may be different in 



various theories of gravity and would help us to distinguish these theories. 

Dynamical evolution of scalar field in Horava black hole spacetimes is studied using 
Horowitz-Hubeny approach [541]. The spectrum of entropy/area is discussed from the view- 



point of QNMs of scalar field in HL gravity 55|, |56) . QNMs of various fields around KS black 



57H60I] . In earlier work 
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we have stud- 
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63|. A 



hole spacetime were calculated using WKB method 

ied time evolution of scalar field using the time domain integration method 
noticeable deviation of the evolution behavior of scalar field in the ringdown region from the 
standard Schwarzschild black hole was found. However the late time evolution of field fails 
to show any distinction from the Schwarzschild case. The late time behavior of massless 
scalar field depends only on the multipole order I and decays as t~^ 2i+3 \ Where as massive 
field in the intermediate range the decays as t~( £+3 / 2 ) sin(mt) , but in the asymptotic late 
time the decay is dominated by t~ 5 ^sin(mt) tail. 

It is interesting to see how various other field perturbations decay in KS black hole 
spacetime. In this paper we study the evolution of electromagnetic and massless Dirac fields 
in the KS black hole spacetime and probe the signature of the now theory by comparing the 
present results with those of the standard theory of GR. The paper is organized as follows. 
In SecJTT] we derive the radial wave equation for various field perturbations around KS black 
hole and discuss the numerical method used. In Sec JIHI the evolution of electromagnetic and 
massless Dirac field are studied using the method of time domain integration and the QNMs 
are evaluated and the results are discussed in SecJIVl 



II. PERTURBATION OF FIELDS AROUND KS BLACK HOLE 

The IR vacuum of pure HL gravity is found to be anti-de Sitter 0, Q]. Even though HL 
gravity could recover GR in IR at the action level for a particular value of the parameter 
A = 1, there found a significant difference between these black hole solutions and the usual 
Schwarzschild AdS solution. The asymptotic fall-off of the metric for these black hole so- 
lutions is much slower than that of usual Schwarzschild AdS black holes in GR. Meanwhile 
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Kehagias and Sfetsos[26j could find a black hole solution in asymptotically flat Minkowski 
spacetimes by applying deformation in HL theory by adding a term proportional to the 
Ricci scalar of three-geometry, /j^R^ while the cosmological constant A\y — > 0. This will 
not alter the UV properties of the theory but it does the IR ones leading to Minkowski 
vacuum analogous to Schwarzschild spacetime in GR. KS solution is given as 26] , 



ds 2 = 

where f(r) has the form, 

N(r) 2 

The event horizons are at, 



N(r) 2 dt 2 + f(r)' L dr 2 + r 2 dVt 



(1) 



fir) 



2(r 2 - 2Mr + a) 



- 2 + 2a + Vr 4 + 8aMr 



r± = M± VM 2 - 



a. 



(2) 



(3) 



When a = the solution reduces to the Schwarzschild spacetime case. 

The evolution of massless scalar field $, electromagnetic field and massless Dirac 
field ip in the spacetime g^ u , specified by Eq.([T|) are governed by the Klein- Gordon pjj, 



Maxwell's 



651 ] and the Dirac equations 66] respectively. 



'-9 



A(v^<TU)$ = o, 



(4) 



F.tr = o, 



with 



flU 



A — A 



(5) 



[ 7 X(^ + r„)] 4> = o. 



(6) 



The radial part of the above perturbation equations can be decoupled from their angular 
parts and reduced to form, 



/ q2 q1 \ 



(7) 



where r* is the tortoise coordinate defined by dr* = jdr and the effective potentials V(r) 
for different fields are given by, 



Vs(r) 



f(r) 



+ 1) Idfjr) 



r dr 



0,1,2,..., 



(8) 
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Vem{t) = f(r) 



+ 1) 



1,2,3, 



V D± (r) 



\k\ — 1, 2, 3, 



(9) 



(10) 



Here k is positive or negative nonzero integer related to the total orbital angular momen- 
tum hy I = \k + \ \ — \ and Vd+ and Vr>- are the super symmetric partners and give same 
spectra. So we choose Vd+ by omitting the subscript. The perturbation equations (Eq.((7j)) 
can be written in the null coordinates u = t — r* and v = t + r* as, 

d 2 



-4 



dudv 



\ff(u,v) = V(u,v)iff(u,v). 



(11) 



In order to get the time evolution picture of the field we integrate this equation numerically 
using the following finite difference scheme suggested in 6^] , 



N 



h 2 

^ w + ^ E -^s- —V(S)(V W + Vjs) + 0(h 4 



(12) 



{V-Vcf 



where the points N, S, E and W form a null rectangle with an overall grid scale factor of h 
having relative positions as N(u + h,v+h), W(u + h, v), E(u,v + h) and S(u, v). Since we are 
interested in the late-time behavior of the wave function which is found to be independent of 
the initial shape, we set ip(u, v — 0) =0 and a Gaussian profile *&(u — 0, v ) = exp 
In all our calculations we set the initial Gaussian with width a = 3 centered at 
order to progress on the above numerical scheme one has to get the value of the potential 
at r(r*) = r((v — u)/2) at each step. Here we use the Runge-Kutta method to numerically 
integrate the equation for the tortoise coordinate and find the values of r{r*) at each step 



2^ 

10. In 



by cubic spline interpolation as suggested in Ref. 

values of ^ on line u = v — 2c is extracted and plotted as a function of t 
the value of r* at which the field is evaluated. 



63| . After the integration is completed the 

v — c; here, c is 



III. TIME EVOLUTION OF FIELDS AROUND KS BLACK HOLE 



The evolution of scalar field around KS black hole was studied in Ref. [61J. The typical 
nature of time evolution of a massless scalar field is shown in FigHJ The plot shows how the 
generic behavior of wave function (initial outburst, quasi-normal oscillations, and power-law 
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decay) in HL theory is differed from pure Schwarzschild spacetime. The QNM phase ends 
at a later time in HL theory and the oscillation frequency and the damping time have a 
higher values in HL theory. We can calculate the QNM frequencies from the numerically 
integrated data by nonlinear \ 2 fitting. The calculated values of QNMs are given in Tabl dlHl 
and TablaHH 
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FIG. 1: Time evolution of massless scalar field around KS black hole with a = 0.8(top curve) and 
the Schwarzschild black hole(bottom curve) for £ = 1. 







WKB 


Time domain 


a 


Re(u) 


Im(uj) 


Rb(uj) 


Im{u) 





0.29111 


-0.09800 


0.29224 


-0.09758 


0.2 


0.29564 


-0.09424 


0.29750 


-0.09418 


0.4 


0.30071 


-0.09008 


0.30154 


-0.09027 


0.5 


0.30345 


-0.08777 


0.30474 


-0.08817 


0.6 


0.30636 


-0.08526 


0.30850 


-0.08564 


0.8 


0.31260 


-0.07937 


0.31416 


-0.07979 


1 


0.31918 


-0.07162 


0.32057 


-0.07247 



TABLE I: QNM frequencies of massless scalar field for I = 1, calculated using WKB and numerical 
integration data 
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a 


Re{uj) 


WKB 

Im{uS) 


Time domain 
Re{oj) Im(u) 





0.48321 


-0.09681 






0.2 


0.49085 


-0.09351 


0.49087 


-0.09356 


0.4 


0.49942 


-0.08969 


0.49911 


-0.08981 


0.5 


0.50412 


-0.08752 


0.50266 


-0.08767 


0.6 


0.50915 


-0.08511 


0.50831 


-0.08527 


0.8 


0.52039 


-0.07926 


0.52029 


-0.07947 


1 


0.53356 


-0.07108 


0.53272 


-0.07138 



TABLE II: QNM frequencies of massless scalar field for 1 = 2, calculated using WKB and numerical 
integration data 



In order to study the evolution of electromagnetic field in HL black hole spacetime we 
numerically integrate the perturbation equation Eq.f llll) with effective potential Eq.(j9]) using 
the time domain method. In Figj2]^a) the time evolution of the electromagnetic field ^(i, r) at 
a fixed radius r* = 10 for £ = 1 is plotted in comparison with the corresponding Schwarzschild 
case. We can see that the time length of the QNM phase increases in HL theory. The late 
time tail starts at t ~ 300 for KS black hole with a = 0.8 whereas it is at £ ~ 230 for the 
Schwarzschild case. The oscillation frequency and the damping time have a higher values in 
HL theory. The variation of oscillatory region of ^f?(t, r) with the parameter a is shown in 
Figj2^b). The oscillation frequency increases with a and shows a slower damping for higher 
values of a. 

The late-time behavior of the field around black hole spacetime is also an important topic 
of study in black hole physics. It was found that during the late-time a massless neutral 
field dies off as an inverse power of time by the factor t^ 2e+3 ^ depending on the multipole 



order of the perturbation [53j, |64J . It interesting to see how the field at late time behaves in 
the HL theory. 

Figj3](a) shows the late time behavior of wave function for different values of a, with 
multipole index £ — 1. We find that the late time behavior is independent of a and follows 
the behavior of the Schwarzschild case with ^ ~ t . In Figj3](b) field evolution for different 
multipole index is shown with a = 0.5. The perturbation dies off at late time as t~( 2£+3 ) 



(a) 




ln(l) 



(1)) 



A/ 



-AAA 

. Ik, 



FIG. 2: Evolution of electromagnetic field for £ = 1. (a) Around KS black hole with a = 0.8(top 
curve) and the Schwarzschild black hole(bottom curve). (6)QNM region of the time evolution for 
different values of a. Curves from bottom to top is for a = 0, 0.4, 0.8 and 1. 

as in the case of Schwarzschild case. The field falls off as \I/ ~ £-5.08^-7.09 an j £-9.09 f or 
£ = 1, 2 and 3 respectively. The predicted values are -5, -7 and -9 respectively. 



(a) 





(b) 



ta(t) 



FIG. 3: Late time decay of electromagnetic field, (a) for different values of a with £ = 1. The Field 
decay as an inverse power of time with t -5 08 , for all values of a. (ft)Decay of field with different 
multipole order £ for a = 0.5. The field decay as t -5,08 , t -7,09 and i -9 - 09 for £ = 1,2 and 3. 



We have seen that quasinormal ringing phase is a dominated form of decay in the evolution 
of perturbations after the initial transient phase. The QNMs calculated from the numerically 
integrated data are given in Tab- Mill and Tabl dHIl We can find a good agreement with the 
earlier results obtained using WKB method in Ref. 60(]. Results show that the oscillation 
frequency and the damping time increase with a. 
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ot 


Re(uj) 


WKB 


Time domain 
Re{oj) Im(u) 


n 

\j 


24587 


-0 0931 1 




-U.UjUtJl 


0.2 


25152 


-0 08942 


0.25225 


-0.08706 


0.4 


25796 


-0 08530 


0.25733 


-0.08447 


0.5 


0.26151 


-0.08299 


0.26383 


-0.08135 


0.6 


0.26532 


-0.08047 


0.26575 


-0.08009 


0.8 


0.27375 


-0.07435 


0.27771 


-0.07454 


1 


0.28308 


-0.06569 


0.28478 


-0.06581 



TABLE III: QNM frequencies of electromagnetic field for £ = 1, calculated using WKB and nu- 
merical integration data 



a 


Re(u) 


WKB 

Im(oj) 


Time domain 
Re{oj) Im(u) 





0.45713 


-0.09506 


0.45088 


-0.09471 


0.2 


0.46531 


-0.09179 


0.46239 


-0.09117 


0.4 


0.47453 


-0.08799 


0.47153 


-0.08711 


0.5 


0.47961 


-0.08582 


0.47835 


-0.08502 


0.6 


0.48507 


-0.08339 


0.48332 


-0.08183 


0.8 


0.49736 


-0.07743 


0.49756 


-0.07439 


1 


0.51197 


-0.06887 


0.52823 


-0.06897 



TABLE IV: QNM frequencies of electromagnetic field for 1 = 2, calculated using WKB and nu- 
merical integration data 

Now we study the evolution of Dirac field in HL black hole spacetime. We numerically 
integrate the perturbation equation Eq.( !TT|) with the effective potential Eq.( !T0|) using the 
time domain method. FigJUa) displays the time evolution of the Dirac wave function ^(i, r) 
at a fixed radius r* = 10, for k = 2. The deviation of generic time dependence of wave 
function in the HL theory from pure Schwarzschild spacetime is clear in the plot. QNM 
region lasts for a longer time in HL theory. The late time tail starts at t ~ 310 for a = 0.8 
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whereas it is at t ~ 240 for the Schwarzschild case. Also the oscillation frequency and the 
damping time have a higher values in HL theory. The variation of oscillatory region of 
ty(t,r) with the parameter a is shown in FigfJJb). The oscillation frequency increases with 
a and shows a slower damping for higher values of a. 



(a) 



A 



I 



ta(t) 




fs. 



4 Art 



ta(t) 



FIG. 4: Time evolution of Dirac field with k = 2. (a) Around KS black hole(top curve) and the 
Schwarzschild black hole(bottom curve). (6)QNM region of the time evolution of Dirac field for 
different values of a. Curves from bottom to top is for a = 0, 0.4, 0.8 and 1. 



FigfS](a) shows the late time behavior wave function for different values of a, with k = 2. 
The late time behavior is independent of a and field decays in the inverse power of time as 
£-5.08 j n pjgj5^b) ; field evolution for different multipole indices are shown with a = 0.5. 
The perturbation dies off at late time as \& ~ £-3.08 ^-5.08 an j ^-7.09 f or ^ — \ 2 and 3 
respectively. The predicted values are -3, -5 and -7 respectively. 

The quasinormal ringing phase of the perturbation can be seen clearly in these figures. 
The calculated values of QNMs from the time domain data are given in Tabl dlHl and Tabl dHIl 
We can find a good agreement of the values obtained in Ref. 59| from the WKB scheme. 
Results show that the oscillation frequency and the damping time increase with a. 

Finally we compare the time evolution of different fields in HL gravity. In Figj6]^a) 
massless scalar, Dirac and electromagnetic fields are plotted for a = 0.5. The only difference 
between evolution of these three fields is in the QNM phase whereas the late time tails follow 
the same decaying pattern with same power law exponent. In Figj6]^b) the variation of real 
and imaginary part of QNM versus a are plotted. All the three fields show the same 
dependence on the HL parameter a. 
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FIG. 5: Late time decay of Dirac field, (a) For different values of a with k = 2. The field decay as 
an inverse power law of time with t -5 - 08 for all values of a. (fr)Decay of Dirac field for different k 
with a = 0.5. The field decay as £-3-08^-5.08 an( j £-7.09 f or k = 1,2 and 3 respectively. 







WKB 


Time domain 


a 


Re(oj) 


Im(uj) 


Re(uj) 


Im(uj) 





0.17645 


-0.10011 


0.17830 


-0.10744 


0.2 


0.18187 


-0.09614 


0.17952 


-0.09481 


0.4 


0.18709 


-0.09167 


0.18756 


-0.09312 


0.5 


0.18963 


-0.08916 


0.19040 


-0.09249 


0.6 


0.19212 


-0.08639 


0.19333 


-0.09172 


0.8 


0.19681 


-0.07981 


0.19635 


-0.07572 


1 


0.20053 


-0.07113 


0.200101 


-0.06939 



TABLE V: QNM frequencies of Dirac field for k = 2, calculated using WKB and numerical inte- 
gration data 

IV. SUMMARY 

The evolution of electromagnetic and massless Dirac perturbations in KS black hole in the 
deformed Horava-Lifshitz gravity is studied in this work. Time domain integration method 
is used to obtain the evolution picture. The question raised was whether we can distinguish 
the fundamental nature of the HL theory from a knowledge of the evolution of fields around 
corresponding black holes. Comparing with the Schwarzschild case, we have found that 
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WKB 


Time domain 




Re(u) 


TmA uj) 

J. 1 t It l tv J 


1 ID I UJ j 


/ Til [/ill 


n 

\j 


37863 


-0 09654 


U.o i o^to 


-U.UooOU 


o 2 


38489 


-0 09317 


U.OOIZU 




0.4 




-0 08932 


0.38935 


-0.09014 


0.5 


0.39583 


-0.08712 


0.39126 


-0.08816 


0.6 


0.39998 


-0.08468 


0.39153 


-0.08388 


0.8 


0.40920 


-0.07869 


0.40337 


-0.07758 


1 


0.41978 


-0.07021 


0.41415 


-0.07239 



TABLE VI: QNM frequencies of Dirac field for k = 2, calculated using WKB and numerical 
integration data 



(a) 




0.45 
0.4 
§ 0.35 
0.3 

0.25 



Scalar 

EM 

Dirac 



ft 
0.2 



ft * * 
0.4 0.6 



ft 
0.8 



0.1 
0.09 
"g 0.08 



(b) 



0.4 0.6 





Scalar 




EM 


* 


Dirac 



FIG. 6: (a)Evolution of massless scalar, Dirac and electromagnetic fields around KS black hole 
with a = 0.5 and 1=1. (b) QNMs of different field as a function of a with 1=1. 

the QNM phase has been extended for a longer time in HL theory before the power-law 
tail begins. Also the oscillation frequency and the damping time have higher values in HL 
theory. The variation of the field evolution on the Horava parameter a are also studied. 
In the time domain picture we find that the length of QNM region increases with a. The 
late time decay of field is found to be independent of nature of field and the HL parameter 
a and follows the same power-law tail behavior as in the case of Schwarzschild black hole. 
These results will help us to check the correctness of the theory provided our gravitational 
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detectors can trap these signals with high accuracy. 
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